A dynamical mean-field approximation study of a tight-binding model for Gai_ x Mn x As 
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The magneto-optical properties of the ferromagnetic semiconductor Gai-^Mn^As are studied 
within the dynamical mean-field approximation (DMFA). A material-specific multiband sp 3 tight- 
binding Hamiltonian is employed for the dispersion of the GaAs host. The calculated density of states 
shows an impurity band and a distorted valence band for large and moderate values of magnetic 
coupling, respectively. Upon using the more realistic band structure, the ferromagnetic transition 
temperature is significantly closer to the experimental results than the previous predictions of k ■ p 
models. The optical conductivity shows a Drude-like peak at low frequencies which is suppressed 
by increasing of the magnetic coupling. 

PACS numbers: 75.50Pp, 75.30.Et, 71.10.Hf, 71.27.+a 
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Introduction- The discovery of a new generation of 
dilute magnetic semiconductors (DMS) with large Curie 
temperature^ ha s le d to numerous experimental and 
theoretical studies!^ The ultimate goal is to find a DMS 
with a ferromagnetic transition above room temperature 
and suitable for use in spintronic devices. One of the 
most promising candidates is Gai-^Mn^As due to its 
rather high T c \^l a nd compatibility with current elec- 
tronic applications. Despite this interest, some of the 
most basic questions about Gai_ x Mn^As are still unan- 
swered after more than a decade of its discovery. Among 
them are questions about the nature of its underlying 
magnetic interactions and the role of the impurity band. 

In this work we combine the sp 3 tight-binding Hamil- 
tonian for the zinc blende GaAs with the self consistent 
dynamical mean-field approximation (DMFA) to study 
the magnetic and transport properties of Gai_ x Mn x As. 
We assume that the manganese (Mn) doping does not al- 
ter the tight-binding parameters of GaAs and any change 
in the band structure is due to the many-body effects in- 
duced by the magnetic coupling between the Mn local 
moments and the itinerant holes. 

Several previous studies have used the tight-binding 
approximation for DMS. Tang and Flatte calculated the 
local density of states for a single Mn and two nearby 
Mn impurities^ and the magnetic circular dichroisnP 
of bulk Gai-zMn^As using large supercells; Masek et 
al. took the Weiss mean-field approach to calculate the 
electronic structure for several DMS^ and Turek et al. 
compared both tight-binding approaches^. Large-scale 
Monte-Carlo studies of real-space tight-binding Hamil- 
tonians including only the three valence bands has also 
been performed^. 

Our DMFA calculation incorporates quantum self- 
energy corrections which were not included in most of the 
studies described above. We also incorporate valence and 
conduction bands on an equal footing while Yildirim et 
ali^ included only valence bands. Because this method 
is non-perturbative, it allows us to study both the metal- 
lic and impurity-band regimes, as well as small and large 



couplings. Since the strength of the coupling between the 
magnetic ions and the charge carriers is comparable to 
the Fermi energy, temporal fluctuations included in our 
DMFA approach are required in a realistic calculation 
especially in the vicinity of the critical temperature. 

We choose ti ght-b inding parameters according to 
Chadi and CoherP^^, which give the correct band fea- 
tures within a sufficiently wide window of ± 2.0 eV 
around the Fermi energy, such as effective masses for va- 
lence and conduction bands and the gap at the center 
of the Brillouin zone. Our model also includes the spin- 
orbit interaction of the parent material which was proved 
very important in previous studies of Gai_ a; Mn :c As!^"^' 
The inclusion of the realistic band structure of the par- 
ent material leads to more realistic results. For example, 
the band repulsion between the conduction band and the 
impurity band, which was absent in some previous calcu- 
lations based on the k-p approaciP™^, results in a signif- 
icant reduction of the bandwidth of the impurity band. 
Narrowing of the impurity bandwidth increases the lo- 
calization of the mediating holes, which suppresses the 
estimated critical temperature. Consequently, we find 
that the calculated T c is significantly smaller, and closer 
to the experimental results. 

Model- We model the magnetic interactions in 
Gai-zMn-rAs using the simplified Hamiltonian proposed 
by Zarand and JankoP^l 



H = H + J c J2S(Ri)-3(Ri), 



(1) 



where Hq is a 16x16 matrix including both electronic 

dispersion and spin-orbit coupling of the sp 3 holes of the 

parent material. Our tight-binding parameters reproduce 

the correct band structure for the heavy and light bands 

around the center of the Brillouin zone, the split-off en- 

er£ 

duction bands 1 

with a term A Q L • S where Ac a and Aa s are interaction 

constants for Ga and As atoms respectively. The second 

term in Eq. (nl describes the interaction between Mn 



rgy gap and the g ap between the valence and the con- 
nnz-.+i™ Konrijn i- The spin-orbit coupling is modeled 
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FIG. 1: (color online). Total density of states versus energy 
at the edge of the valence band for x = 0.05, nj, = x/2 = 
0.025 and several coupling strengths at T—58 K. The edge 
of the conduction band lies at the right bottom corner. The 
impurity band is not well separated from the valence band up 
to large couplings ( J c ~ 4 eV). The vertical dashed line shows 
the location of the chemical potential for J c =4.0 eV. Band 
repulsion between the conduction and impurity bands leads to 
reduction of the impurity bandwidth. Inset: sp 3 total density 
of states for pure GaAs and Gao.95Mno.05As with J c =4.0 eV. 
We use a Fourier transform filter to reduce the noise in the 
calculated density of states. The noise is in regions away from 
the Fermi energy and is due to the effect of the k mesh in the 
coarse graining step of the DMFA self-consistency loop. We 
have used up to 62,000 k points in this calculation. 
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FIG. 2: (color online). Average hole polarization versus tem- 
perature for different values of J c at nj, = x/2 = 0.025. The 
polarization is frustrated since it does not reach its maximum 
allowed value, V = 3/2, at T — 0. Comparison with Figure 
[l] shows that noticeable changes in T c happens when the im- 
purity band is forming. Inset: Average magnetization of the 
ions versus temperature for the same couplings. 
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spins and mediating holes, where J c is the exchange cou- 
pling between the localized moments with spin S(Ri) and 
the holes total angular momentum density J(i2j), both 
at site Ri. Since the spins of Mn +2 's are relatively large 
(S — 5/2) we treat them as classical vectors. 

We employ the DMFA algorithm^ to calculate the 
magnetic and transport prop erties of the material within 
a Green's function formalis m 1 17 ! 18 !. We calculate the av- 
erage magnetization of the manganese ions: 



M = 1 J <m s S z exp{-S eS (S)} 



(2) 



with the partition function Z — J dQ, s exp[-S cS (S)], and 
the effective action S Q g(S) = — J2 n logdet[GQ l (iuj n ) + 
J C S ■ J] J22B2 where Go is the cluster excluded Green func- 
tion. 

The average polarization of the charge carriers at the 
chemical potential is defined as: 



V 






(3) 



where jf is the z component of the total angular momen- 
tum and rii, the density of holes in the ith band. 

To compute the optical conductivity we use the Kubo 
formula: 22 



where Tr is the trace, Uj(k) is the velocity along the i 
direction, A(\l,uj), the spectral function at momentum 
k and frequency ui, and /(w), the Fermi function. The 
group velocities involve the hopping terms in the Hamil- 

tonian and are calculated as Uj(k) - 



SA, 



A=0 



where 



H is the Hamiltonian and A the vector potential. 

Results and Discussion- We first focus on the den- 
sity of states (DOS) of the charge carriers at the doping 
x=0.05 where T c is among the highest reported! 1 ! 23 ! 24 ! In 
Figure [T] and [2] we assume 50% hole compensation due to 
anti-site and interstitial doping, such that n^ = x/2 = 
0.025. This corresponds to the optimum filling in the 
impurity band regime. In any case, the profile of the 
DOS changes minimally with filling. Figure [T] shows the 
total density of states at the edge of the valence band 
for several values of the coupling. Most experimental 
probes, such as photoemissiorpS, infrarecpSHSU and reso- 
nant tunneling 2 ^ spectroscopies, and magneto-transport 
experiments^ infer a value of J c between 0.6-1.5 eV. 
Within this range of couplings the effect of the magnetic 
interactions is just a distortion at the edge of the va- 
lence band. We have to use coupling strengths as large 
as 2.5 eV to observe the formation of the impurity band. 
However, recent scanning tunneling microscopy experi- 
mentsSH display an impurity band similar to the one ap- 
pearing in our results for couplings around J c ~ 3.0 eV. 
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FIG. 3: (color online). T c versus filling at constant doping, 
x~5%, for different couplings. For small couplings T c de- 
creases by increasing the hole concentration. For larger cou- 
plings there is an optimum filling of n^=2.5% where T c is 
maximum. For coupling values lager than 2 eV, T c does not 
change appreciably. The inset shows T c vs. coupling strength 
for rih = x/2. The saturation of T c is due to the lack of 
non-local correlations in the DMFA. 
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FIG. 4: (color online). Optical conductivity in the direction 
perpendicular to the magnetization for different couplings at 
x—5%, nh— 2.5% and T—58 K. The conductivity is suppressed 
by increasing of J c due to the increasing localization of the 
holes. The inset displays the total density of states versus en- 
ergy at the edge of the valence band for the same parameters. 



The impurity band is not completely separated until cou- 
plings around ~ 4 eV. As we expected, the repulsion be- 
tween the conduction and the impurity band confines the 
latter within the band gap and narrows down the impu- 
rity bandwidth. This leads to localization of the charge 
carriers which in return reduces the ferromagnetic tran- 
sition temperature, T c . 

Figure [2] displays the temperature dependence of the 
hole polarization for different couplings. We choose a 
narrow energy window (~0.1 eV) at the chemical poten- 
tial and calculate the average V within this window using 
Eq. ([3|. The polarization is frustrated since it does not 
reach its maximum allowed value, V = 3/2, at T = 0, 
and it is further reduced with decreasing magnetic cou- 
pling strength. Comparison with Figure [T] shows that 
noticeable changes in T c correspond to couplings where 
the impurity band is forming and separating from the 
valence band. On the other hand, T c does not change 
once the impurity band is fully separated from the va- 
lence band. The maximum T c we obtained is ^220 K 
which is closer to experimental result s^ than previous 
estimates using the k ■ p modelPZES One can also see 
the same behavior from the inset which shows the aver- 
age magnetization of the Mn +2 ions versus temperature. 
However, note that within the DMFA the average mag- 
netization is not frustrated. This is an artifact of DMFA. 
Inclusion of non-local correlations using cluster methods 
reduces the value of the magnetization!^ 

Figure [3] displays T c versus hole filling for x—5% and 
different couplings. For large J c , in the impurity-band 
regime, T c is maximum at the optimum filling of n^ — 
x/2. This is consi stent with studies based on the im- 
purity band picture!- 14 ! 17 ! For moderate couplings around 



J c ~1.5 eV and fillings n^ > 2%, T c increases by decreas- 
ing the hole concentration. Recent studies in thin films 
find T c to be proportional to the hole concentration for 
a wide range of n^ valuesP^H This might indicate a fill- 
ing larger than the doping, nn > X, in those experiments. 
The inset in Figure [3] shows T c versus J c for fixed hole 
concentration nh — x/2 — 2.5%. The plateau at large 
couplings is due to the absence of non-local correlations 
within the DMFA. 

Finally, we calculate the optical conductivity in the 
direction perpendicular to the magnetization, er xx , for 
energies smaller than the band gap and different cou- 
plings (Fig. EL For small couplings, the conductivity 
displays a Drude like peak which is suppressed by in- 
creasing the coupling. This suppression is due to the 
increasing bounding between the magnetic ions and the 
itinerant holes. Our results are in agreement with other 
theoretical calculations. 9 However, we do not capture the 
low energy (~0.2 eV) peak observed in measurements of 
the conductivitypO The origin of this peak is still highly 
controversial. While Burch et al. explain it as evidence 
of the existence of the impurity bancF§J, others explain 
the result within the distorted valence band picture.^! 
Since our results for the density of states do not show 
any low energy feature for J c <2.0 eV, we do not ex- 
pect to see any peak in optical conductivity for J c <2.0 
eV. Moreover, for large couplings the wide impurity band 
will induce a feature with similar width. The inclusion of 
non-local correlations reduces the bandwidth of the im- 
purity band and, it might be possible to capture this low 
energy peak within a dynamical cluster studyP^l 

Conclusions- We study the magnetic and trans- 
port properties of the diluted magnetic semiconductor 



Gai_ x Mn x As within the framework of the multi-orbital 
DMFA. We employ a semi-empirical sp 3 tight-binding 
approximation to model the band structure of the parent 
material. We choose tight-binding parameters according 
to Chadip^l which reproduce the correct band structure 
within the relevant energy window, including the valence, 
split-off and conduction bands. The spin-orbit coupling 
is modeled with a term A Q L ■ S. We find that this more 
realistic band structure leads to more realistic results. 

The density of state shows a distorted valence band for 
moderate coupling strengths, J c ^1.5 eV, and an impu- 
rity band for couplings larger than J c ^3.0 eV. The band 
repulsion between the impurity and conduction bands 
reduces the bandwidth of the impurity band. This in 
turn leads to localization of the itinerant holes and re- 
duction of T c . The average hole polarization displays a 
clear magnetic frustration. Moreover, the optical con- 
ductivity er xx shows a Drude-like peak at low frequencies 
which disappears with increasing couplings. While our 



conductivity results are similar to previous tight-binding 
calculations,^ they do not display the low energy peak 
observed experimentally!^ Despite the inclusion of band 
repulsion effects and that our impurity band is narrower 
than found previously, it is still too wide. This wide 
impurity band smears out any low energy features. We 
believe that the inclusion of corrections including cavity 
field and non-local correlations into the dynamical mean- 
field approximation will result in further suppression of 
the impurity bandwidth and, possibly, development of 
the peak in the conductivity. 
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